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Abstract. A study of the effect of deformation and pairing on the development of halo nuclei 
is presented. Exploratory three-body core + n + n calculations show that both the NN interaction 
and the deformation/excitation of the core hinder the formation of the halo. Preliminary self- 
consistent mean-field calculations are used to search for regions in the nuclear chart where halos 
could potentially develop. These are also briefly discussed. 


MOTIVATION 

During the early years of Radioactive Nuclear Beam physics, while the dripline for light 
nuclei was being explored, nuclear halos became a very hot topic [ 1J. Typical examples 
of halo nuclei include 6 He, u Li, 14 Be and 19 C, all in the low mass region of the nuclear 
chart. The big open question, when moving toward heavier systems, is whether halos 
can develop when for instance A > 40. 

In order to successfully describe a halo nucleus, the structure model needs to take into 
account Jj]: i) the very low density region in which the halo nucleons move, subject to an 
interaction that is closer to the free NN interaction than the realistic in-medium nuclear 
interaction; ii) the long tails of the wavefunctions and correct asymptotics of these tails, 
which contribute decisively to many nuclear properties; iii) the few-body dynamics of 
the valence nucleons relative to the core and between themselves. For these reasons, it 
is acceptable to decouple the halo degrees of freedom from the core’s, simplifying the 
standard microscopic treatment: this is the basis for applying few-body models to light 
nuclear halos. 

In heavier systems, the decoupling of core and valence nucleons may not be as 
straightforward. Mean-field studies in the past decades have shown that pairing is im¬ 
portant, and in some regions of the nuclear chart (namely the deformed region), defor¬ 
mation is also necessary. As matter densities in Hartree Fock Bogolyubov (HFB) have a 
faster radial decay when pairing is included, one expects the pairing force to act against 
the halo formation. The pairing anti-halo effect was shown in Ref. |:5] through self- 
consistent Hartree Fock Bogolyubov calculations for Carbon isotopes. In opposition to 
[5], the results of Hamamoto et al. [6,@] indicate that when an sq / 2 neutron approaches 
the dripline, it becomes decoupled from the mean field, which means they are less paired 
and consequently may give rise to a halo. One can then claim that, if the neutron becomes 
decoupled from the mean field, the HFB theory should not be used. Unfortunately, the 
calculations in are not self-consistent; the mean field is reduced arbitrarily, from a 
well bounded situation, to force the system toward the dripline. 


If the situation with pairing is controversial, the effect of deformation appears to 
be well settled. Nilsson model calculations performed by Hamamoto [j|] suggest that 
J z — l/2 + states (J : being the projection of the angular momentum on the deformation 
axis), which contain s-wave components of the valence particle as well as other compo¬ 
nents with higher partial waves, become pure s-waves when the system it forced toward 
threshold. Thus, it is concluded that deformation does not hinder the halo formation [11. 

This exotic feature is not solely associated to nuclear physics: halo manifestations ap¬ 
pear also in atomic, molecular or condensed matter physics. A detailed review [2] com¬ 
piles the results across the nuclear border, and develops a halo signature that removes 
the specific scale (nuclear, atomic, etc). 


RESULTS WITH THREE-BODY MODEL 

In Ref. fl] we study the structure of the ground state of a nucleus with two valence 
nucleons as the system approaches the two particle threshold. A three-body model of 
core + n + n is used where the core is deformed and allowed to excite. The Faddeev 
Equations are solved within the hyperspherical method. Our starting point is the 12 Be 
model [§]. An important ingredient of the model ]§] is the core+n interaction fitted to 
reproduce the properties of 11 Be. We artificially decrease this core — n interaction, 

Vn-coreCr) = A V^Jr) . (1) 

with A —> 0 to simulate the approach to the neutron dripline |4J. 

The correspondence from a few-body language to the mean-field terminology appro¬ 
priate for heavy systems, is far from trivial. Nevertheless, in some way, pairing and 
deformation are included in this three-body model. In the limit of a very heavy core, the 
core-n interaction plays the role of the mean field in the microscopic description. The 
NN interaction V nn included in the three-body model would then be related to the pairing 
associated with the valence pair only, in the mean-field language. However the pairing 
force and V nn are not identical. Here V nn is a sum of gaussians fitted to the low energy 
NN phases shifts, containing central s, p and d terms, as well as a spin orbit and a tensor 
force, while pairing in a typical mean-field calculation is zero range and s-wave only 
(see 0, IlOO for a further discussion). 

The hyperspherical expansion introduces a new quantum number K related to the 
sum of the angular momenta associated with the Jacobi coordinates (x,y) |9]: K — 
(4 + /y)/2. The resulting coupled channel equations contain a centrifugal barrier of the 
form (AT + 3/2)(^-|-5/2)/p 2 , where p = \Jx~ +y 2 is the size of the halo. Even for the 
lowest hyper momentum K = 0, there is a barrier reducing the halo effect. For a halo to 
appear, the K=0 component needs to be the dominant component (>50%) [2j. Thus, we 
look at the probability of K=0 in the ground state wavefunction as a function of the three- 
body binding energy (Fig. 1). The probabilities for the lowest hyperspherical harmonic 
component (K=0) as a function of the two-neutron separation energy are shown for four 
different cases: i) the case where a realistic NN interaction is included, as well as a 
quadrupole deformation = 0.67 for the 10 Be core that reproduces the experimental 
B(E 2;0 + —> 2 + ); ii) the case where only the NN interaction is switched off; iii) the case 
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FIGURE 1. Probabilities for the K=0 component in the ground state wavefunction of a 12 Be-like system 
as a function of the two valence neutron separation energy. 


where deformation/excitation is neglected but the NN interaction is included; iv) the 
case where both the NN interaction and core deformation/excitation are neglected. 

The first thing that should be noted is that, even for no deformation/excitation and 
for no NN interaction, the system never develops a pure K=0 state in the limit of zero 
binding. This result is associated with the fact that we are treating the system as a three- 
body system and it is in contrast with what would be obtained in the two body case 
0]. Secondly, both the NN interaction and collective core degrees of freedom decrease 
considerably the limiting value for Prob(K=0) when Sin —► 0. Ultimately, when both 
effects are included, the probability is well below 50 %, suggesting that no halo will 
appear. 

Note that the physical input for this exploratory study was provided by 12 Be. One can 
repeat the calculations reducing the deformation parameter and increasing the mass of 
the core, to simulate a heavier system more realistically [if]. However the model would 
still be inconsistent, in the sense that the core-n interaction would not be determined 
from the physical subsystem, and the approach to the dripline would remain artificial. 


PRELIMINARY HFB RESULTS 

The appropriate way to address the problem for medium mass to heavy nuclei is through 
a self-consistent model that contains both pairing and deformation 0. Let us first 
concentrate on pairing only. A systematic HFB study was performed, searching for 
nuclei with the 3.s ( y 2 
0. All calculations 

found that the Cr isotopes satisfied this condition. Using SLy4 IUJ for the particle-hole 
channel and ULB zero-range density-dependent pairing [[l5f] for the particle-particle 
channel, the neutron dripline is reached for 80 Cr. We plot the proton and neutron Helm 


neutron orbital intercepting the Fermi level in the limit of stability 
were performed with the HFB code written in a 3D mesh IH .We 




radii [16] and geometrical radii (Fig. 2) as a function of isotopic mass. The large 
difference between the neutron and proton Helm radii indicate a neutron skin. A neutron 
halo exists when there is a large difference between the neutron Helm and geometrical 
radii. We see that both are present on the Cr dripline. A detailed analysis of our results 
do not confirm those from Ref. |6] concerning the decoupling of the s-wave orbital from 
the core, when reaching the dripline iH. 


Helm criteria 



FIGURE 2. Helm criteria for halos and skins: geometrical and Helm radii for both neutrons and protons 
in Cr isotopes. 


Helm criteria, z axis 



FIGURE 3. Same as Fig.(2) for Dy isotopes along the deformation axis. 

Similar searches where performed in the region of deformed nuclei using HFB with 
even multipole deformation IH . For the deformed region the orbital we need to concen¬ 
trate on is the harmonic oscillator l/2 + Nilsson orbital coming from the 4 s |/ 2 shell in 
the limit of sphericity. Only for the prolate Dy isotopes did we find this valence neutron 
orbital to be close to the Fermi level at the dripline. The analysis based on the Helm and 














geometrical radii was performed on the Dy isotopes (see Fig. 3), similarly to the Cr case. 
For that purpose, an extension of the Helm criterion to deformed systems was proposed 
0. As there is no significant difference between the neutron’s Helm and geometrical 
radii, we expect no halo to develop. Note that Dy isotopes demonstrate another minimum 
for oblate shape, which means that the correct ground state should have configuration 
admixture. Such a calculation was not performed due to its complexity, yet we do not 
expect it would reverse the result concerning halos. 

From these preliminary self-consistent HFB studies we can draw the following con¬ 
clusions: 

• In a realistic situation, pairing does indeed oppose the development of the halo. 
This confirms the results of o] an d the three-body model [4]. The effect is not 
sufficiently strong to hinder completely a halo in 80 Cr. 

• The self-consistent HFB results do not show the decoupling of the neutron from the 
mean field seen by 0]. 

• Deformation in heavy nuclei seem to strongly hinder the appearance of a halo, 
as suggested by the three-body model [0] and in disagreement with [8]. No halos 
were found in the deformed region of the nuclear chart. In particular, the fully self- 
consistent calculations do not confirm the purification of the J z — l/2 + orbitals at 
the Fermi level into s-waves when approaching the dripline. 

An article containing all details of these HFB studies will become available soon. As part 
of our near future plans, we will study the dependence on the choice of the microscopic 
pairing force, namely a better low-density behaviour and finite-range effect 0. 

• This work has been partiaffy supported by Nationai Superconducting Cyciotron Laboratory 
at Michigan State University and the National Science Foundation under grants PHY-0456656 
and PHY-0456903. 
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